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ABSTRACT 


A commonly studied problem in the field of cryptography 
is (i Discrete Logarithm Problem. This problem is also 
referred to as the “distance” problem. Basically, one would 
like’ to know where a particular piweryaseusie 2S: 1 a list 
combining all of Neti, Bepecs cuted as ‘powers of some 
primitive ST atienty. or Saud eaicntiy what is the distance - 
between a given mr of on-tuples in a \similar 
representation. A de Bruijn sequence is 4: well-known . 
| periodic binary sequence in which every n-tuple from 0 Es 
2°-1 appears. Our goal is to beErey understand es oe 
ones” de Bruijn eceusnee. ite ae eh to understand 
“where each of the binary n-tuples eee in that panei. 
Using the Necklace Algorithm, the Sequence of n-tuples can 
_ be generated. This list has some Speer ee operates that 
allow us to perform the required analysis to locate the 
n-tuples by an association io. nesses: We partition ‘the 
binary aeuoke ie necklace classes accérding te the 
“longest substring of ones appearing on the n-tuple. We nen 
| count how many n-tuples appear ies ope Sequence for the first 
time as members of a necklace class containing HG aienGer 


strings of ones. 
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I. INTRODUCTION 


A. BINARY N-TUPLES 


» The Ser .0f binary n- tuples is the set of all strings of 


“IOS arid l’s of length nm. Usually hese are ordered from 


O00. “he: 11s These ae eupites can be -oune, on necklaces 


defined by cyclic rotation. | These Hemeiaecs are Lyndon 
words of length d, where din. Example for nm = 4: 
(0 ' (0001) (0011) (01) (0001) (1) | | 
These are the Lyndon words (0) and (1) of length one, (01) 
of length two, and (0001), (0011) and (0001) of length four. 
B. NECKLACE ALGORITHM — 
Pe. Definition of Necklace Algorithm 


The Necklace Algorithm is an algorithm for generating 


the necklaces of beads of length n in two colors, i.e., the 
binary necklaces. Fundamental to the Necklace Algorithm is 


; the operation, 0: y OV pas defined as: 


© (ar, Aoypeup An) = (bi, b2,..., Dn) 7 


where b; = a;.for i = 1,2,..,j-1 and j is determined as ‘the 


largest subscript such that a>0 and a, = 0 for ail. k>j. Then. 


ob = a1, and by: = b, for ft =1,2,..,n-j. 


Necklace Algorithm (for necklaces of length n). 


Step 0. The initial necklace, V,', is 11.1 = 1"; 






Step 1. OQ@-step: Form ov”) COIN Vi a7 





Step 2.  J-check: The resulting string ‘is the next 






necklace if and only if jin. 





Step: 3. If we have not found the last necklace, 0”, 





by steps 1 and 2 above, return to step 1. [1] 





Figure 1. Necklace Algorithm 
Note: If step 2 is omitted, we produce all Sener Laces 
(defined fetenh. If step 2 is modified such that mite we 
produce a list of all Lyndon words of iengen on. 
The following output of the Necklace ‘icowtvn ange 
the three options for the J-check is provided for n ae 


No J-check 





Table 1. Necklace Algorithm Output (n = 4) 
Thus the Necklace Algorithm produces a listing of the 
n-tuples appearing from the largest to the smallest and with 


dof the n-tuples appearing on each Lyndon word of length d. 


2 








PAdteeonetig, “fhe: wenden words that form the Sceiascs en 
‘be juxtaposed in oe: to create the de Bruijn “prefer- 7 
ete? sequence as follows: £5111 021100.104100020< [2] ue 
2.7 Properties of the Necklace Algorithm 
The Necidtaee Algorithm euepue is the list. of n-iong 
necklaces, where the largest n-tuple on Sach: necklace 
represents the necklace and all the npee aces appear in’. 
“lexicographic order from largest to smallest. Other authors 
order necklaces: end their respective Lyndon words from 
smallest to largest and represent a ngekiace. Be Sts teese 
8 eepeeen otis. “ie esate are eqgieicns and it is easy 
ey translate between them. Inherent in this list, — Mead 
n- tuples are partitioned into classes defined according to 
the length of ‘the initial string of 1’s in the n-tuple. We 
make some Bierier ve definitions. 7 
Definition 1. An n-tuple of length 7 has ee 
; ae | | eS n is even , | 
parameter I, where [= a | , and is a member of the 
i 7 eee is odd | : 
class, m. - These parameters also indicate the number of 
leading 1's as there are I-m ieading ivs pErone the first 0. 
Each class then has the Fneyat form with the “respective 


parameters as follows: 





om. DOF | : 
1...10@,_,.5+.-@, such that l-n < m < lI. 


_Note: The extreme special cases are the n-tuple, 1”, 
where: m=l-n since there is no zero following the lead 


string of ones, and 0”, where m=I1 since no ones preceding 


the first zero. 
| l—m 


..The leading ee is defined ae Peng Oiy’ and is : 
common to all members of a respective m class for a given l/l. 
To consider the number of neubeee ie a class, we also 
consider the degrees of Passe (DOF), i.e., those bits > 
that are unspecified by the characteristics of the ease: 
We consider a particular n-tuple class as an eeempies 


V = 1110s 670g j 0011 


[: _ 5 

n (Odd): 11 4 2h1 = 2(5)+#1 = 11) 

m: 2 (e.g. I-(I-m) = 5-3 = oF 

DOF: 7 (e.g. 2l41-[ (I-m)+1] = 11-[ (5-2) ye) = 7) 


A subtle point is to be noted. The leading 1’s in a 
class m will increase in number with /, and thus with 2, 


while m remains fixed. So an 11-tuple in class 3 has two 


leading 1’s, while a 13-tuple in class 3 has three leading 




















“1's.” The il- ~tuple has 8 degrees of freedom while the 
13- tuple has 9 degrees of freedom. | | | = 
The general form of the eco, of the Neektiece 
Algorithm is refined ‘to represent the oucout ‘into . the 
respective m etacees: ‘The first OUtPUE.'1n class m srpears 
by applying © to the last Soteik of the class m-1. The 
Sake Sieuenk x of the class m-1 is ee as (i-m+1] ane ‘the. 
first element of the class m is then denoted as © [l-m+1], 
‘where — 7 


Lh net | 


“Jom l-m —di=m n(enod(l—m+1)) 
i meamthomeaan’ 


efi - m+ 1j=11...1011...1011...10... 111... , 
and the last output of the Necklace Algorithm for the class, 


m, referenced as “ats hes ne form: 


l-—m 


114000... a 

i Prenecklace Generation we 
B prenecklace | appearing by step 1 of the Nedieiace | 
pene is defined as follows: a _—— 
‘Definition 2. An aovedtex, y", is a prenecklace 5 Eee 
there exists 2éne Beeccar y* such that yr = Pa ae yy ; eicKe 
i yn is a nevis of length m+k. | (Note: £0) ; denotes the 


operation of concatenation between V" and M@) [3]. 








Theorem 1. The vector V” ee is a peeneciies if 
and only if Wi... Va-Ke1 > VpuVn for all k such that S Isksn. Bae 

The Lyndon words of length d for din, when exeenasd a 
“length n, form all the necklaces of length n, The Lyndon 
words of length d ford < n, when extended to length a Poti | 
all the prenecklaces of length n. [4] We describe the 
relationship between the prenecklaces, necklaces and the 


Lyndon words in the Table 2 and Figure 2 below for n= 4. 


Lyndon Words 





Table 2. Comparison of Type of Words 


choke. ds a weer mapping of the — of all Lyndon 
words of ‘lengths one through n pace oe set of all 
prenecklaces of length n BY taking each of the Lyndon words 
of i d < nm and Ceres them to length n. These see. 


then the outputs, as defined by @, to make the n-long 


prenecklaces. [4] 

















Ruskey shows that to count Am), the — cae 
prenecklaces of Length n, we need only sum a ene 


l<i<n, aig of the se Words of smaller or aGael Length: 


P()=SL@. 7 2203) 


i=] 


Equation 1 
“Likewise there is. a Bev eceren from che: Set of “all 
7 ee words of length d, such that d|n, onto the dhieian 
Of aes n. To count mann the necklaces of Length n, ve 
- need only to sum all of the Ey noe weeds oi length d such 
that d\n: 


¥,)= v1, @. 


din 
epaEOns 
‘These relationships imply that the necklaces are a. 
subset of the prenecklaces. This is illustrated in the 


following figure: 





Figure 2. Neek lace Miger cums Reduction of Bosecs Space 
Prior to the publication of the Necklace Algorithm [1], 
generation of the necklaces was seeouelisnes by selecting an 
eae and performing n cyclic shifts to determine the 
representative ee the necklace class. Graphically, this 
would be similar to picking a point in the cee. and 
moving from point to point n times until one point was 
determined that is the representative of the necklace class, 
N’. By step 0, the Neeiiees Algorithm begins the search EOL 
necklaces in WN’. | By @, he search for the next necklace 
continues in ee set of prenecklaces, P’. > [3] _=&By 
constraining our search to this subset of all n-tuples, we 
eearipe - savings in the uae: tag strings that are 
considered as a possible neciace., Then the J-check removes 
the need to perform the eveiie rotation to test whether the 


n-tuple represents the necklace. 








Teableats. demonstrates oe ‘alternative ' methods a 
counting nee eeneciaeee of length nm. . First the abd ed tie 
based on the number of Lyndon Oras of length i, 7 bli) ; is 
es [6] Uesnani TS, the eeisigice Algorithm partitions te 
prenecklaces into their respective m classes. Table | 3 is 


provided for the case n = 11. 


it 

w ‘ 
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Table 3. ‘Comparison of Counting Prenecklaces 


A Seite prenecklace will appear in different 





locations in the two tables. (e.g., 1%! is located in the 
rows i = ] and m = -6, while 0*' appears in row i = 1 and 
m=5.) The discussion of how to count the prenecklaces 


with respect to the class m appears later and is the major | 
contribution of this thesis. By carefully counting the 
nae in the Necklace Algorithm caliente we determine ‘the 


location of any pap hece ee n-tuple in the nrerere ones de 





Bruijn sequence. The savings neseure for each m class and 
its difference from the power of 2 described by the | 
respective degrees of freedom is the number of aeeiaecs 
appearing as enumerated ai tha second part of Table 3. 

4. Necklaces | 

By step 2, the J-—check eeinoues those prenecklaces that . 
are not necklaces. From our statements in the previous 
section (section 3); step 2 therefore removes the Lyndon 
words of length jin. | This results in additional savings to 
those observed for step 1. The following ereorss were ies 
the claim fiat ss 2 produces only necklaces: | | 

Theorem 2: A vector, VW", satisfying step 1 of the 
Necklace Algorithm, | isa neeiEes if and only if j|n. 

Proor. ° ee ae vector, = —— euepose n= fq. 


then vy V2.uV_ = (Vie 0)'. Clearly this is a necklace only 


when y, --V,40 > Vv, Vig 


Ov,...¥,, Wi¥0- 
Now let mn=f+s with Is<s<j, and assume without loss of 
Ov, ...Y;) Wie0, as before. 


generality that y,...v,,0>y,...v,, 


Then y,...V,¥,...V,,0>¥,...V,_,0Y...V,7 since dropping the first s | 


bits in each string yields Vz oY 0 > Vey oY), 0,...¥, and the 


prenecklace is not a necklace.p 
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. Alternative to Equation 2 above, we also. have 
MacMahon’s formula for the number N(n,2)of binary necklaces: 


NG, 2)=— Ly gaan : | | | | 17] | 


Pain 


Equation 3 

Here Ad) is Euler’s totient function snd Cher seen 2S 
over all divisors of ehetaseieee length, n. 

ks mentioned in the discussion of prenecklaces above, 
_ the formula does not partition the necklaces into classes as 
does the Necklace Algorithm. | tenes: the Necklace 
Algorithm does not enunceaes the necklaces in classes as we 
i would like. (Note: see Appendix A) 
a dt an analysis to compute the size 
- of the classes for this parameter. From the definition of | 
sud tasece: eye only ¢ one necklace class nae a 1 An 
the last position, namely sae Any other prenecklace that. 


ends ina 1 can simply be evel teetty rotated to the right to 


apes the 1 in the first position: making the resulting 
vector aie than the panes vector. e.g., 1101 = a EG, 


and 1110 >/1101. 
From the previous See teiente: an elenentary Eppes bound 
for the number of. members i. ae class can be determined. 7 


Since the iast “Position of a necklace must be a QO, an upper- 


ie 








bound is 2°. the savings described above shows how far | 
below this power of 2 the number of necklaces actually is. 
By careful analysis of the Necklace Algorithm, this upper 


bound can be tightened. 
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sds COUNTING THE NECKLACE ALGORITHM 


A. PRENECKLACE COUNT 
1. Table of Data for Missed n-tuples 
Our ultimate goal is to better understand the “prefer- 

ones” de Bruijn sequence. We wish to show where each of the 
| biases eet pies appears in the sequence. The first step is 
to Nebaiins to which necklace the n-tuple belongs. For 
ssch necklace it ie eeeetene to determine its class, as the 
binary n-tuples have been partitioned into si asdeatuecoruins 
to the longest Neen substring of ones ‘by the Necklace 
Algorithm. If we can count how neny necklaces appear in the 
- sequence for the first time as members of a necklace class 
containing he, tenes strings of ones we will nis eouna te | 
location o€ the n-tuple. We choose, in fact, to count nen 
‘many necklaces fail to appear dia, eee way (i.e., those 
necklaces that appear earlier in the sequence as members of 
other necklace classes containing longer strings of ones as 
well). That is, our anna sie proceeds as we count cher 
- strings nae are ‘missing “from the m-class ee. Same aides 
iets: see Table 3 and subsequent discussion) . 
The computer code nee to implement the Necklace 3 


— Algorithm is modified to count the number Of n-tuples 


i3. 





omitted (in the above sense) from inclusion . in the 
respective m class in the list of prenecklaces. dee. he 
respective m classes the number of possible n-tuples with 
n-(l-m-1) degrees of freedom is re NO aul of these 
appear employing the Necklace Algorithm as we have noted. 
The number of mere is used as a performance measure for 

the efficiency of the algorithm as stated above and in 
Ruskey. [3] The entire table of data ais included in 
Appendix A. These tables in Rppendis A, as presented, are 
pegnegetes by odd and even values for the iensth of the 
n-tuples. Our analysis will be for the case in which nm is 
odd. When eC even the sutivetc is similar. The data in 
Table 4 and Table 5 is erevecesd from Appendix A. This data 
specifies the nee of n-tuples missing from the list of 


prenecklaces of length m in the class m. 








fal 7[s [ups 









[ea [== [253[_ 489] 20] 704 
res [== [== [2023[ 1930) 2877 
Tes [a= pa | == [2098 7526 
0 ae ee ae ee ELE 


Table 4. Number of Missing n-tuples (Odd) 
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sel 58 
Ces [= [== [anil sae aass 
5 We al ad LL LL 
per pepe pt par 


Table 5. Number of Missing n- SupESs (Even) 




















Table 4 ane Table 5 provide several insights EOI our 
further analysis. First, there are no missing n-tuples when 
the number of leading ones is greater than or equal to l, 


that is, when m < 1. Second, when there are ‘no leading 


ones, i.e., m = 1, the only appearing n-tuple is 00...0° 





Thus, there are 2" 2-1 n-tuples missing that have the err 
string of one’s being Ses. honey For example for n 
L=3, the last n-tuple from the Necklace Algorithm belonging 

ie the re m = 3 is 000000, which — I-m=1-3=3-3=0, 


i.€., the leading number of ones is zero bai ie degrees 


of freedom ae is. five. and 31 n-tuples are omitted in 


this class. The dashed lines in the tables indicate that no 


value is defined for ae given parameters as the numbers are 


not meaningful. _A third observation is that the numbers in 


the table seem to increase to a value and remain at that 


ees 








value. We call this the steady state number for that m 
etaas. The values when steady state has been reached are 
indicated in bold. This meeeeay seats number reflects the 
maximum savings afforded by performing 4n6 Necklace ~ 
Algorithm for that class. faeiiedeiy, the frest Suh eeoien 
we face is to count the size of Cae wtendy state set ae 


missing n-tuples for a given class m. Further, we show for 


which value of hn the steady state occurs for a given class 


m. 

The notion of a class -reaching steady state er 
demonstrated in Table 6. Geastaee the following n-tuples 
from the class m = 2, which are missing from the list of 


prenecklaces for the lengths of n: Ty 9, and ll. We do not 
include the value of = 5, which for the class m=2 has 
zero leading 1’s. Thus there are 15 missing members in the 


class, namely the non-zero 5S-tuples with a leading 0. 


16 














n=7 n=9 | n=i1i1 Count 

1011111 = 110111111 LITO TTI TI i 
TOO os OTTO 11101111110 2. 
1011101 = ~° 110111101 11101111101 3 
1011100 -s 110111100 LITO 7100: * 4 
1011011 s 110111011 © 11201111011: 5 
oar. => 110110111 11101101111 . 6 
1011010 = 110111010 11101111010 q 
1011001 xs 110111001 137021717001 8 
1011000 = 110111000 11101111000 Q 
OLOldtT. as “12010112 11101011111 10 
1010110 s 110101110 11101011110 11 
-1010011. => 1170100111 : 11101001111 12 
1001111 s 110011111 11100111111 .. 13 
1001110 s 110011110 11000111110. 14 
1001101 s 110011101 11100111101 15 
| =  . 4110011011 “2E1 00111011. ., .. *-16 
1001100 = #£110011100 11100111100 © -17 
1001011 s 110010111 11100101713 18° 
1001010 = °»=110011010 11100111010 19 
1000111 = 110001111 11100011111 20 
7000110 x= © 110001110 11100011110 24 
- 1000101 x=. #£4«~110001101 11100011101 22 
1000011 s 110000111 11100001111 23 

21 23 23 | 





Table 6. Increasing Missing n-tuples to Steady State 
oat n= 7 the m = 2 class e not yet ss steady state; 
n=9 ie the first value of n that the class m a7 is. “t 
steady States At n= 11, tee Necklace Algorithm produces a 
fist: of necklaces that remains at steady state and there are 
| twenty-three n-tuples omitted for the class m= 2 for each 
eyes cain value of n. Comparing the iceiel: leita. + whe dhe -— 
that each eaine aeene can be cnr (=>) es me teeta” 
in ene: successive column by adding a 1 to the iéeaing 


substring (this is required LOX to maintain the class 
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structure), and adding an additional 1 to the longest string 
Gf I's in the trailing sees Ge eee ee n= 7 to 
n = Q, there occasionally are two opportunities to add the 
second 1... It is interesting to note that ‘from n= 9 to 
n=11 each missing n-tuple maps to one and only one ‘eats 
an the meet column. For enampie, 1011011 produces > Boch: 
F101 011. ana: JOT 01d. i wanton 3 we soe whe such 
possibilities for two successors end and steady epece is 
therefore achieved. | | 

2% Plot of Missing n-tuples Relative to m 

An approximation eG the Peon number of missing 
n-tuples is made by analyzing the data in Appendix A. We 
extract from the data the first tag th n for each class m 
at which steady sees is saeaiel oe analysis includes the 
values for the parameter n (2741 for odd and 2/ poe euch l, 
m, i-m, the degrees of freedom and ie, eae BunbeE of 
missed n-tuples. | The aanbee of ‘n-tuples missed is 
transformed using the logarithm base two to aid in the 
scaling on the graph. (Base 2 seems a natural enous Since 


we are concerned with a binary alphabet. ) 
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21] 1 | m | I-m | dor | #missed | log0 | 
i ee ee Re ee” eee: See: 
| 9 4f 2] tT 28] 4524 


rasa 4)_io] 138] 7.109 
ai, tof al | _i4]___740[__ 9.5311 

Peis s-__s|__ie 3720] 11.861) 
[33,16 of dof 22) 17936[ 14.133) 


Table 7. aEeacy _ for odd values of n | 





aT a [or ateses Toa 


7 : 12 - : : 
1672[ 10.707 
sof is] 7) 82a] _8208| 13.003] 
36 ie] e| 10 25| 38944] 15.249 


Table 8. Steady state for even values of n 





ieee ae ag identiy a detest ieee between land 
the selene of missing n-tuples. The following seapnec Ok i 
vs. log2(#missing n-tuples) illustrates this nearly fee 
relationship. For reference, a line of the form y = x is 


plotted in dashed lines as well. 
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Odd Steady States 


a, 
a=) 
® 
2] 
ae 
E 
Ra 
— 
N 
ae | 


Log2(#missed) 





Figure 4. Regression of Steady States (Even) 
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This linearity is iegieaidl: adi a reer regression of 
the data where the seeeeted sum of Sauaucs yields _" 
correlation of determination (R’) of 0.99790 for ‘the odd 
length n-tuples and. 0.9989 for the even Sep ies. this 
ends ua to propose 2! ~ steady seaee value. However, it ro 
also evident that enis is only a loose upper bound as I (and . 
therefore Ay eueue: our further analysis will also serve to 
tighten this eee | 

3. Finding Steady State 

‘The eens: Algorithm proceeded using m as a 
parameter. § The above argument establishes / as a 
determining factor for the number of missing n-tuples (which 
is also pe oa linear relationship with HH). FOr the odd values 
of aA, eeeey state is reached for on = 3,9,15,21,27,33;«. 
Evidently steady state is reached when n=3 (modulo 6). : We 
ee structure of the n-tuples when the Steady state 
is Peaened: The predecessor and the first n-tuple in the 
class where steady state is reached exhibit a similar 


property among the classes. 
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a v2 
ne ©(111000000) rn 110110110 
pbb lo ile Sibiemoneds J =] fnlonnor1110 | 
ae ee de 
pe fao far @(f11111100000000000000) m/l firitionti01n1 110 


Table 9. Examples of First Member of Steady State Class 






By the Necklace Algorithm, ite age one in an n-tuple 
is complemented producing I-m-1 ones and a zero. Le, 2s 
apparent that the-length, Hm, of this sinensis poekers 
divides nm with a cofactor Of 34 “We also consider: larger 
values for n es ensure that steady Stace: as maintained. 
Table 10 “illustrates the first sutput of che, Necklace | 


Algorithm for the class m = 2 for various lengths of n. 


re 
~-@(1100000) | = #10101 a 
i-] 2 
aie . j 
ST @(41i000000) 10110110 
| 2 
pa os ioooon | | 11101110111 
pas | os ~ (4117700000000) | {111011110111 


Table 10. i of First Member of Class m = 


















. 











it can be sane iad coon the data in Table 4, for the 
class [~2, steady state is reached at 2H] = 9, and continues: 
for all heeaey ‘vata of ‘, The first keine of the : 
Necklace Algor thn for the class m «a aaa fon eae value 
2 of n or waiew steady pees is reached nes a following - 


form ce is used to abbreviate the remainder) : 


1-2-2 R- 
| emt TET meetin) 


1...101---10ll@, . 
7 The term remainder iat cates ‘that Part of the vecaine 
substring chat is not completely i laa | Evidently, this 
shows elise if f is large enough ee. allow the leading 


substring to be ‘copied two © times and leave a remainder oS, 


| length three or greater then eee. state will Be reached a | 


for the gies m= 2. The eiecsee of larger values Be m | 
also exhibit a distinct relationship with 1 and m. . For 


ne 21+1= 31, and the classes of m = 1, 2, 16 the imei 


similar neeaetenieties apes Ras Table 11). 
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TSE S| 
class [- class /]-2 Remainder 
eer 













14 


—— : 
100000) —i 1., jo cer 


cas 

Paya 100000) = dic 101. f...do1111111 

i—6 1-6 ) 
oxi. io0000) | >| 


py maa 
1...101...101...11 
Table 11. Examples of First Output for Classes (J = 15) 























(Note: for m = 6, steady state has not been reached for 
Rie= Si.) | | | : 


This sugeeses ene eeition ae eeoren: 

Theorem 3. 

Steady state for the @-step of The Newiees Algorithm 
will be cnaciend for the class of m at on = 2H1 ’ when 
1+22>3m. | _ | 

Proof: From Table 10, we see that he sunsering “110” | 
is copied three times for 2l i l= 9. For larger n 4, ee 


[=m 


same m, the initial string 1°-"0 will only be copied twice 
plus the Heese 2m-1 bits This establishes the following 


inequality 
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3(/-m+l)>2+1=n 
31-3m+3221+1=n 
[+223m> 
To illustrate Theorem 3, we verify the result for the 


ntuple of length 31, showing that the class of m = 6 and 


the number of leading ones (l-m = 19-6): is not at steady 


state: 





Table 12. Verification of Theorem 3 
(Note: Steady state is achieved at 2I+1 = a8) 

By the theorem, steady state is suesinabic, "Because 
hie is not a priori apparent Ge: oce this specifically in 
the following eoeslieeies: | | 

'» Corollary 1: 
For a given class, m, che Necklace Algorithm will peach 
steady state for some I > a. 
Proof: m a6 a fixed quantity. Suepece that. ae, 


Let L = 3m. Increase J such that 1 > L. QO. 
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Corollary 2: 
‘Seed State is reached for spuds G2 and only pe ee 
yan ie with a sailed: of 2m-1. 
It is helpful to consider _ example se Corollary 2. 
Let / = 15, then nv = a1: From Table LZy we can compare the 
results for the classes m = Lee She 54.0 asin the 


corollary. 





(Note: Steady state is reached for all cases but m = 6.) 


B. NECKLACE COUNT 


£. Table of Data for Omitted Necklaces and n- Piedes 

For the “savings” realized by the a we fivet 
examine the difference between the prenecklaces and ive 
neereees: of length an. | eaaaeien E and Equation 2 
respectively rae values of P and N’, respectively. 


These values are presented in Table 14. 
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Table 14. Prenecklaces, Necklaces, and the Difference .. 
The number P” is aponexmaeeny twice the suite: a 
more information can be obtained by comparing the number of 
prenecklaces and necklaces for each si ite deat aaa for a 
fied Wi: In particular, steady see becomes apparent > 
again. Table 15 sieee values savtitioned by ‘m-class ocr | 


the case nm = 31. 
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im Fi) | Nin) Nim) | Pin)-Nim Nia 












oj a2 768] —ressal “Tessa 

—esss3[_s276e| a7 
ay ast0aat—ess22} 65827 
[6 2079218] 1038850] 1040368 
[sf 8005876] 3990960] 4014916 
-—9]15076085| 7492056] 7584029 
a CE | | 


Table 15. Count of Step Outputs (n = 31) 




























Tt apparent that P(m) -N(m) grows by a factor of 2 until. 
m > 1. Then when the Necklace Algorithm is poesonrrct a 
necklace representative is found without error the no 


rotations and comparisons. But the cost of the algorithm is 


still about twice the number of necklaces foun Since the 
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eC eee have to be. determined first. peeiewine Figure 
2, the number BE eecieei ices: that are not necklaces, 
P(n)-Niy), he eoieeneie on the order Ss the number of 
| necklaces [ais Thus the work to determine the necklaces is 
cut dois ‘from the n2". rotations and comparisons’ to 


OW) =02"/) as a function of utilizing the Lyndon words. 


C. DEFINING STEADY STATE 
Steady state has been shown to “Occur for particular 
m-classes of Sidnucktunes in reaceds: 3. In | the previous 
witht va prenecklaces are determined to ee n-long 
extensions of the Lyndon words of length d, for 1 < d <n. 
Similarly, the necklaces are determined by extending to . 
length ye Lyndon words of length d, where din. tie new. 
consider partitioning the inden words die m-classes os 
verify thes nerien of steady state. | 
From the partitioning into m-classes of the 
- prenecklaces for a given ny, _ can ‘count che number of 
sEGHSe chaGes . ier a ee length of n using Equation 2. | 
‘Since the prenecklaces are d-long Lyndon words extended to 
iength nm, we can count the number of prenecklaces in a ae 


class by counting the number of d-long Lyndon words that 
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oe with the substring 10; Thus we can enumerate the 
Lyndon words for each m-class and for any value of nv. 


Table 16 lists the number of prenecklaces of length nine 


determined by Lyndon words of length d< 9 and partitioned | 2 


by the number J-m of leading is 





Table 16. Partitioning Lyndon Words 


(Note: The ney eee ee for l-m represent the 
ieeecGa aun order aiposed by the Necklace Algorithn. The 
Necklace oe proceeds by columns, from left to acne | 
in the Rabie: Bacau? each n~- ead extended LyneoP word 
(prenecklace). | 

From the table ae is evident ‘that the number of 
prenecklaces in the class ee is (141424+3+6+10418 = ae. 
In Tani. 4: Seeady state for the class m = 2 is 23 missing 


n-tuples We verify that | 2bOF -Pin) = 64 - 41 = 23 The 
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reader can easily check that steady state is reached for the 
class . = 1 (i.e, missing 3 Caples in the column | 
iy rn en | 

~ Not only can ‘it be determined noM any prenecklaces are 
missing for a eee anee class, but it can be determined how 
many are missing for a particular value of d. | By , 
" precalculating the size of the set be aesing. n- -tuples at 
| stecay efeee for each m-class, it ais possible. ‘to precisely 
‘determine the position for a given neeriaee. A necklace in 
a class Lor which steady state has been achieved can easily 
be calculated using the known value for a precise position. 
_ For ‘other classes, not at steady state, this information can 
“Be epee in eotapieenind tight upperbound for phe? dumber 
of necklaces in a given class before inieiaeing a sequential 


search in.the list of necklaces. 
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- III. CONCLUSIONS AND SUGGESTIONS FOR ADDITIONAL RESEARCH . 
_ The Necklace Aigosithn- “nae “been shown ~ to ; be 
Somieiceney more efficient ‘een al naive algorithm 
Bpeaueine: all of ene necklaces of a given Ronen a 
teviosanaenie order. The Seen cnn also has been shown to 
*» treaties all of the Lyndon oe of iehethe one ecusn oe 
lexicographic order as prenecklaces. — By developing ine 
notion of steady state, it is shown to be possible to find 
the mack position of the necklace representative of any 
particular n-tuple given that the acing substring has at 
least I ones. — pe: the representative Sualave has fewer - 
leading 1's, a table is available to count the missing 
necklaces for the given class m. The class ere ce at 
| steady state to ensure that the pabievsize e feasible. The . 
condition for steady State is shown to _ be 1+2<3m. 
Additional methods must be employed when steady state as 
not een sen ieued: | | — 
To. demonstrate the application of our results, - - | 
| example is given. We choose an n-tuple: _ 
pp = 000110111110101. 
First, we determine the n-tuple’s parameters using 


Definition'1. The longest substring of 1’s is 11111; and 


og 
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‘since 1 = 7 we find m = 2. Now, we check for the steady 
state condition: 


742 > 3(2) > 9 > 6 > Class m = 2 is at steady state. 


From Appendix A, the steady state values for odd values > 


of m are 3, 23, 138 in the classes m = 1, 2, and 3. The 


representative of the necklace class is: -111110101000110. 


This representative has no internal periodicity so the > 


J-check would have shown that this n-tuple isa Lyndon word 


of length d= 15. There are 242 Lyndon words of length 15 


and with l-m = 5. - To locate “the exact necklace then we 
would enter the string [i-m-1] ( = 111111000000000) . - We 
then apply the Neekiaee, Higceitan and count how peers 
necklaces are produced, until the correct one is éouaa. 

| From this example and our discussion in the previous 
section, | it is clear that a better Ghaerse nada of the 


Lyndon words and the development. of an efficient counting 


technique for a partitioning by class would result in an 


improved solution for this version of the distance problem. 


With the knowledge of the location of a necklace ‘relative to 


an initial position in the “prefer- ones” de Bruijn sequence 


a distance can be determined [5] and refined quickly Ln such 


a manner that it would be computationally efficient. [8] 


This could result in the development of a new cryptographic 


34 











scheme. 
in data 


LOe Unis 


Our methods and results may also have application 
compression. Certainly, it can be said the solution 


problem is valuable for both academic and practical 


reasons. 
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APPENDIX A. TABLE OF MISSING N-TUPLES COUNTS | 













(3 | 63] 104] 128[ 136] 138] 138] 138] 138] 138] 138] 138] 138) 138 
4 | [25s] 459[ 628] 704] 730|_—738|_~— 740] 740] 740) 740] 740] 740 
5 | |__| 1023] 1930[2871{ 3392{ 3606] 3684]_—3710| 3718] 3720) 3720/3720 
6 | | | [409s] 7926] 12584] 15631] 17038] 17606] 17822] 17900] 17926] 17934) 
7 | ft Tf 16383) 32157] 53717] 69872] __78079| 81646] 83100] —-83670|__—_—83886) 
jes [Tt [Tf 65535] 129653] 225253] 305619] 349970|_370379| 379118] 382732 
ro [ [ tT | tt [262143] 520955] 932995] 1315402] 1542563] 1652306] 1701131) 
mio {| Tf tt tf 1048575] 2089242| 3830876] 5592057] __ 6710766] 7275055, 
uf of 4194303) 8369683] 15630420] 23543815) 28890516) 
m2 {| of fT fT tT 16777215}33508823| 63477929] 98359882) 
mis [| [ | of tt ft 67108863) 134107119] 256902720 
nia [| | fT tT tf 268435455] 536601228) 
PS Se ad oe | 


Table 17. Missing n-tuples (n odd) 
















("T7729 [a [as as 7 a7 [a9 7 21 7 23 2s 27 2s 
zo | | of of ol} lof 





4 | | | 204] i69] 76] 26] 8] CT 

ws [| | ff 907] oat s2if a4} 78] 26] 8] tT 

ne | | tf 3831) 4658] 3047] 1407] 568] 216) 78] 6] 

a7 ft ft fT 15774] 21560] 16155] 8207] 3567] 1454] 570] 216 

re tT | ft tT 64118] 95600] 80366] 44351} 20409) 8739] 3614 
|| 1040667] 1741634] 1761181] 1118709] 564289) 
tf 4175380] 7260737 7913395} 5346701) 
Lt *f6 731608] 29969106} 34881953] 
po 66998256) 122795601) 
po 268165773} 
ee re Fae eee See CS 


Table 18. Forward Difference (n odd) 
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| te Ts fof 12] a4 | a6 | 18 | 20 | 22 | 24 | 26 | 2 | 30 
vo | of of of of of =6o) ll} 
5 CS | CO 7 | | (epee | S| ee | Se | | | ee | | 
ee ee ee ee a) ee ee ee) ee ee | ee) ee ee) 

13 | 31) 48] 56] 58] 58] 8] ss] 8] 58] 5) SS] SB] 58 
4 [ [27] 221) 288] 314] 322] 324) 52a] 324] 524] 324] 324] 524 
es [ [| 51a} 946] 1353] 1558] 1636] 1662] 1670| 672] 1672] 1672] 1672) 
| tf 2047] 3920] 6038] 7319] 7878] 8094] 8172] 8198] 8206] 8208 
wf ot CT 8191] 15983] 26075] 33160] 36599] 38044 38614) 38830] 38908] 
rw oT fT ft | 32767] 64609] 110209] 146483] 165716{ 174327] 177932|_179388 
zo [| [ ft |) ft ff 131071] 259975] 459035] 635182] 736127] 783764804579! 
vio | | | | | 524287] 1043452] 1892322] 2715891] 3222158) 3472217) 
Ht ft tT ft ft 2097151] 4182089] 7743350] 11486697] 13940496 
fi2{e[ [ | TT ft | 8388607]16747871] 31514779] 48164396] 
His {| ot | UE 88554431] 67037905] 127749326 
4 | of ft 234217727] 268262880) 
es ae Ea ae a | eS ee 


Table 19. Missing n-tuples (mn even) 






























es {| | | 435] 407] 205} 78] 26] 8 eae) 
re | | of tf 1873] 2178] 1281] 559] 216] 78] 6] 
7 of tT tT 7792} 10092] 7085] 3439] 1445] 570) 216] 78) 
re [ | ft ft ff 31842} 45600] 36274] 19233]__861i] 3605] 1456) 
ro | TT Tt tT 128904) 199060] 176147| 100945] 47637|_ 20815 
Hio [ [ | Tf ttf 519165] 848870] 823569] 506267) 250059 
fir [ [ [. jf tt ff 2084938) 3561261) © 3743347) 2453799) 
iz [ { | tf | 8359264) 14766908] 16649617 
m3 [tT | fT tf) 38483474) 60711421) 
Aa re = 5 nero se oe eae oe eed MO, 

NG Me 2a ce — Ue i ll en cee ee ca aie ene ed 


Table 20. Forward Difference (n even) 
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APPENDIX B. CODE FOR NECKLACE ALGORITHM 
“The following code is a Java pEegren that implements 


the Necklace Algorithm: 


import java.util.*; 

import java.lang.*; 

import java.io.*; 

public class Theta { . 

- public static void main (String args[]) { 

int length = 31; 
int. k= 1; | 
int sequence [] = new int[length] ; 
int j = length-1; 


int counter = 0; 
int numNeck = 0; 
int numType = 0; 


int flag = 1; | 
int group [] = new int[lengtht+1]; 
int sum = 0; 
System. out.println("Finding all binary necklaces of Reweeh my 

length + "."); 
// Sets the necklace of all 1's. 
for (int a = 0; a < length; at++) { 

sequence [a]=k; 

System. out. print (sequence [a] ) 
} 
System.out. printiln(* Whe 
// Find last non zero element and decrement it 
while (sequence[0] != (0) { 

ij = length-1; 

while (sequence[4]. == 0){ 

. J=)= 1; 
_ // Change last one-bit and fill register. 
- seguence[j] = sequence[j]-1; 

if (j !'= length-1) { 

for (int 1 = counter+tl; i < length-j; itt) | 

sequence[jti] = sequence[i-1]; 

} _ | | | a 
if (length @(j+1) == 0) { | // JT-check 
| for (int a = counter; a < length; at++) { a 
' - System.out.print (sequence[a]); 


} 


3g 
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APPENDIX Ci canter LYNDON WORD COUNTED BY a 


Lm 


a A 39 sf 0 
J OE OO 
A 
yaaa 38a ane 25a aoe] 9a] asa] 22m] ae] 710] 0 
ea] a[ ae] eae] 255] Sal 909] asa 3180] anes] 2083] 36h 0 
[ee 3[ ea aaa 286] sia] Toro] — inea| seta] enzo] 7a] size] 92] 0 
[a0] te] 32] ee] aa] 356] sin] Tore] —a0i3[ ~soi0] —7ase| —vatia| — 13239] — soca] 7300 
[| 33[ 64) anal a5 — sn] 1003] aise sore] 776] —vasea| —asete| —aaeaal 1028] Tica 9 
[2a [ 64] Tae] 256] i] 1025] a] — aoeil — s00e| isasal —zasse| — aso75] —sevei|“24563|_T7I| 0 
a A a A Bl Be Ba Ws We Bi in 
[a4 [986] S12] Tons] 007] 050] 8157] —ieise| —s1889[_ei07a| —Tr0ssa] “Fras07| isa7as| w209| A305] 0 
[as | sia|_vezs|aosa|— dona] —ai7a| —res0el sx5s2| Geo0] iaieza| nese] —sao0sa] _see3s6) 15e59e| 67709 
[on aoa in Tease aaa eT 
saat] e120] as02a| Bsz4ea| 7900] _¥50524| Tas07IG| TSTRIse| SOLS 16264] 0 
Ec a 
Ee tM Le ar a aki ane Balen Bo ln 
30 | 1s] 52762] ess01| Te] Bora] siDR00| oa7s5azonLesaaTeReAs| eseRaT7| Sse7eRe| WSGESI [28977] 61567] 0 
a a a 


| Table 21. n= 31 (part I) 
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ae 
| a a 
EC 
ee ee 
i a ee 
Fe OD ee 
BR BR a sR a 
so [fof af af a 48] te] 52] 64] aaa] 256] 51] i024] 2048] 4056 
rap oop ap ayaa ae aa al al ea Toe Soe ane ae] Tee 


Table 22 n = 31 (part ITI) 
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